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Abstract 

We present the calculation of the next-to-leading corrections to the tree amplitudes which appear 
in the description of non-leptonic B-decays in the factorization approach. These corrections, together 
with radiative corrections to the jet functions, represent the full next-to-leading contributions to the 
dominant hard spectator scattering term generated by operators Oi,2 in the decay amplitudes. Using 
obtained analytical results we estimate B — > nn branchings fractions in the physical ( or BBNS) 
factorization scheme. We have also found that the imaginary part generated in the hard spectator 
scattering term is rather large compared to the imaginary part of the vertex contribution. 
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Introduction 



The processes of non-leptonic decays of B mesons are considered as one of the most interesting topics 
at present time. They are sensitive to the physics of the standard model and provide a nice possibility to 
search for new physics effects. The channel of two 7r-mesons is very interesting because phenomenological 
analysis of the corresponding branchings and CP asymmetries can be done to a good accuracy in a 
model independent way pQ. Last few years B — > ttt: branching fractions and CP-asymmetries have been 
measured by BABAR [2] and BELLE [3] collaborations. 

Important progress has also been achieved in the theory. There was suggested a new approach which 
is based on the idea of QCD factorization. The factorization allows, in some sense, to constrain the 
strong interaction background in a model independent way and therefore provides a theoretical basis for 
analysis of B— decays which can be considered as an alternative to the traditional phenomenological fits. 

The factorization theorem for nonleponic decays has been initially suggested in [4]. The statement 
has been proved by explicit calculations at the leading and next-to-leading orders. The general proof of 
the factorization to all orders can be done using the so-called soft-collinear effective theory (SCET) [5]. 
The application of SCET technique to the two mesons decays has been formulated in [6] [7] . The presence 
of two hard scales /x ~ mf, and /i ~ yAm^ leads to two steps matching QCD^ SCETi — * SCETn with 
corresponding two independent coefficient functions which can be calculated systematically in perturba- 
tive QCD. The non-perturbative dynamics is described by the matrix elements of the light-cone operators 
constructed from the fields of the SCETn effective theory. These unknown functions are universal for 
all processes and therefore can be constrained from the global analysis. Phenomenological applications 
to the nonleptonic decays based on the factorization have been considered in several papers. The QCD 
factorization approach (the so-called BBNS or physical scheme) was used in [S]. A different analysis on 
basis of SCET was suggested in [8] . Although both approaches are based on the same theoretical idea 
they are different in the consideration of some phenomenological moments, see for instance discussion in 

An important question which appears in application of the factorization is applicability of perturbation 
theory at relatively moderate scale /j, ~ yAm& ~ 1 — 2 GeV . This situation arises at the second step of 
the matching SCETi — > SCETn • In order to answer this question, the next-to-leading calculations of 
the so called jet coefficient functions have been done in [12l [13l [14] . It was demonstrated [13] that the 
radiative corrections are large but, on the other side, do not indicate any problem for the applicability 
of the perturbative expansion. But the full next-to-leading contributions also include corrections to 
the hard coefficient functions which describe matching QCD to SCETi effective theory. A priory, 
such corrections could also be considered as a source of quite large contributions, especially for the 
color suppressed amplitudes in the BBNS analysis [13]. Therefore the second tail of the next-to-leading 
corrections, corresponding to the matching of QCD to SCETi at fi ~ also have to be computed. 
An other important motivation for such a calculation is the observation that the imaginary part of hard 
spectator amplitude arises only from the radiative corrections. If it can produce sizable corrections to 
the CP-asymmetries then such contribution is very important for the phenomenological analysis. 

Recently, such calculations have been carried out and results are presented in [15] for the graphical 
tree amplitudes and for the penguin amplitudes [16 . In this paper we present the calculations of the 
radiative corrections to the graphical tree amplitudes. Our results have been computed using different 
technical approach and can be considered as an independent derivation of the corresponding corrections. 

Our paper has the following structure. In Section I we introduce the basic notation and review, for 
convenience, the formulation of the factorization theorem for B — > rnr decays. In Section II we discuss the 
matching from QCD to SCETi . We define the basic set of SCETi operators and recalculate the leading 
order coefficient functions. The calculation of the one loop diagrams and results for the coefficient 
functions are given in Section III. Section IV is devoted to the numerical estimates of the branching 
fractions. The discussion of some technical questions and the analytical results for the individual diagrams 
are presented in the Appendix. 



2 



1 QCD factorization for B —> ittt decays 

For the convenience we review shortly the basic QCD factorization approach suggested in [9]. The 
amplitudes of two pion decays are given by matrix elements 



A^ = (n(p')w(p)\H eff \B(P)) 



with the effective Hamiltonian 



G 



10 



H eff ^-^\^ \ c l O u l +C 2 0^ + Y J C l O l + C 7l 7l ^C Sa 8g \+h.c 



G 



i=3 
10 



\ CrOl + C 2 O c 2 + C i°i + C Tr°7i + C 8g 8g j + h. 



(1) 

(2) 
(3) 



where a£^ = V u bV* dl — V c bV* d , Cj and 0\ are coefficient functions and effective four-fermion 

operators respectively. In particular, the explicit expressions for the current-current operators are 



Oi = (u b) V -A{du)v-A, 02 = i u fi b a )v-A(d a Uf3)v-A 

Q\ = (c b) v ^ A (dc) v ^A, 0§ = (c/3 b a ) V -A(d a cp) V -A 



(4) 
(5) 



where as usual V — A = 7 M (I — 75) and indices a, (3 stand for the color. The definitions of the remaining 
terms are standard and can be found, for instance, in [17j . Taking into account the structure of the 
effective Hamiltonian © the decay amplitudes can be conveniently rewritten through the effective 

ng way [i 

Ar+»- = -\i d)l -^-M 2 B .U [a, + al + at EW ] - AW^jMliW [a% + o^ tEW ] , 



amplitudes on in the following way [9] 



npl 

i 



(0) 



A^ = xi^ M y 7T 



-a 2 + a 4 - -a 3 EW - -a^ EW 



l 3,EW 



l 4,EW 



(7) 

where we have neglected the annihilation contributions. The amplitudes cnj describe the matrix elements 
of the different operators in ([3]). Namely, 011,2 gives the matrix elements of the current-current operators 



0i,2 ; ol a ' and 



u,c 



denote the QCD and Electro- Weak penguin contributions respectively. From 



the isospin symmetry one has 



V2A n o n - = Afford + A n +n- 



(8) 



We used notation Mb for B— meson mass, / ff is pion decay constant and below fo = fo(0) — /+(0) 
denotes B — ► 7r transition form factors at q 2 = 0: 



(n(p) \qrb\B{P)) = f+{q 2 



P^+p"- 



q 2 



(9) 



The amplitudes include all dynamical information about the decays. In the limit of large 6— quark 
mass TOb — * 00 the QCD factorization approach makes it possible to calculate amplitudes to the 
leading power accuracy. Let us consider the matrix elements a\ t2 which provide dominant contribution 
to the branching fractions. Their expressions are given by 



Oii=fo / du Vi(u)(p„(u) + / du (p„(u) / dz Ti(u,z) ^(z), 
Jo Jo Jo 

/>oc i*l 

^{z) = JbU \ duj I dx (/> B (u)J{z,x,w) <fn(x) 
Jo Jo 



(10) 

(11) 



We have slightly changed the original notation removing fo from the normalization factor 
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where functions Vi and Tj are the hard coefficient functions which can be computed in the perturbation 
theory order by order in QCD coupling asll 

V t = V^iu) + ^Vf LO (u,z,m b /fi F ) + ... (12) 
Ti = Tt°(u) + ^lf LO ( M , z, mt/fip) + ... (13) 

Z7T 

The hard coefficient functions describe the hard subprocess in which quarks and gluons arc highly virtual, 
with typical hard momenta p\ ~ m 2 . Performing integration over such fluctuations we reduce QCD to the 
effective theory SCETi which however still contains large hard-collinear fluctuations of order pf ic ~ m^A. 
Integrating over these degrees of freedom we reduce SCETi to the low energy effective theory SCETn 
which contains only collinear and soft particles with small p 2 ~ p 2 ~ A 2 off-shell momenta. The coefficient 
function which appears at this step is the so-called jet-function J{z,x,u>) : 

J(z,x,w) = a s (fJ.hc) J LO {z,x,u) + ... (14) 



where the hard-collinear scale \ih c ~ \/ m^K. The soft physics is encoded by the matrix elements of 
SCETn operators constructed from the soft and collinear fields. These matrix elements are parametrised 
by non-perturbative light-cone distribution amplitudes (LCDA) tp v ,<pB and decay constants fg. 
Their explicit definitions are given by 

Utp«(x)=i J ^e- < ( te - 1 )^ A <7r-(p)|J(An)fi76«(-An)|0) ) (15) 

where n and n are the light cone vectors: n 2 — n 2 — 0, (n • n) = 2 and pion decay constant defined as 

(tt-(p) |3(0) ^75 u(0)\ 0) - -if v (p.n) , (16) 
that implies J dx ^^{x) = 1. B— meson LCDA is given by 

FstatitiVM^M") = -ij ^ e+luX <° |g(An)^ 75 ^(0)| B{P)) , (17) 

where v — \{n + n) is the velocity of B~ meson at the rest frame. The Me-independent decay constant 
F s tat(n) [II] is defined as 

Fstat(fJ-) = f B /K F (ii), K f (jj) = 1 + ( 31n y - 2 ) ( 18 ) 

where the physical decay constants /b is given by 

(0 |?7 M 75&| B(P)) = if B M B v». (19) 

As it was shown in |20[l21j. the normalization integral for 0b (w) and higher moments are not defined and 
therefore the non-local matrix element (jTTJ) can not be reduced to the local matrix element (fl"8"|) . This 
feature makes this function quite different from the standard LCDA 's of light mesons. 

As we can see from equation (|10[) the jet function appears in the second term only. This term describes 
the hard spectator interaction. For the fist term in (fT0|) the matching to SCETn is not possible due to 
the overlap of the soft and collinear regions, see discussions [THUS]. Such contribution is known as soft- 
overlap form factor. In the BBNS prescription this form factor is excluded using the so-called "physical 
scheme" . In this approach the soft-overlap form factor is rewritten as a sum of physical form factor /o 
and of hard spectator scattering contribution (the details are given below in the text) . 



3 In this paper we always assume that perturbative expansion of any quantity R is defined as R = R LO + ^-R NLO 
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Explicit expressions for the hard coefficient functions read [9] (i ± 1 = i + (— f 



Vi(u)= [a + 



^Vr°(u) + 0(al), 

Z7T 



V; 



NLOi 



N r 



-V(u), 



rn b 



V{u) = -C F 12 In — - f 8 + 3 



In u — iir ) + 



2Li 2 [u] - In 



21nu 



1-w 



1 - 2m 

1 - u 

(3 + 2iix) In u — (u «-> u) 



rp ( \ Cj±l 1 0! Sr r NLOf \ , rtf 2-1 

Ti(u,z) = — + — (u,z) + 0(a s ), 



N c 1-u 2?r 
For the jet function in our notation we have 

CV S(x — z 



J(z,x,uj) = a s 



N c x umib 



^J NLO (z,x,u) 



(20) 
(21) 
(22) 
(23) 

(24) 
(25) 



The next-to-leading order expression for J NLO (z, x, w) has been recently obtained in several papers 
[T2l [T3l H3] and we shall not present it here. 



2 Calculation of the hard coefficient functions 

The aim of this section is to discuss some details relevant for our calculation. We shall reproduce the 
leading order results for hard coefficient functions Tj quoted in eq. ([M)l . As it was discussed in the previous 
section, Tj is associated with matching QCD to the effective theory SCETi . Technical details of such 
calculations have already been discussed in [23l [24J [25] for the case of heavy-to-light currents. 

First, let us fix the basis of SCETi operators relevant for our case. We introduce two operators with 
approprite flavor q and chiral structure: 

J {A °\s) = (q n W c ) s (l - h v , (26) 

y i n d J 

J {B1 \s) = — {q»W c ) s (W^iIp ±c W c )^ s) (l - l5 )h v . (27) 
nib 

where we have accepted the notation introduced in [26) . The light quarks are supposed to be collinear 
fields in SCET approach fiq n (x) = 0, h v is HQET field. Notation W c is used for the hard-collinear Wilson 
line involving only n ■ A c component of the collinear gluon field: 

(q n W c ) s = q n (s n)Pexp \ig I du nA c [(s + u) n] > . (28) 



Matrix elements of these operators between physical particles define two SCETi form factors: 

( tt+(p) \(d n W c ) h v \B{P))= m b (29) 

( 7r+(p) \j( B1 \ s )\B(P)) =m b J dz e <-(2*-i)™» £\ z ), (30) 

where dependence of the form factors on mass m^, is implied. In order to obtain the factorization formula 
(fTO"]) one nas to perform matching of the effective Hamiltonian to the operators in SCETi (|26[) and 
([2"7| . We shall focus our attention on the contributions of the current-current operators Ox^ because 
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they provide dominant part of the two body decay amplitude. Then for the matrix element A vv ([I]) we 
obtain 

Ai„ = ((tttt), IdO? + C 2 0- 2 \ B) = ~-^M%U{ a % ), (31) 

oii = ( £tt «i * <p x + tp„ * *i * £f *) , (32) 

where Uj and U denote hard coefficient functions and by asterisk * we denote, for simplicity, convolution 
integrals. The index i is introduced to distinguish two possible final states (7r7r) i=1 = 7r + 7r _ , (7T7r) i=2 = 
7r°7r°. Corresponding matrix elements define amplitudes a\ and a 2 respectively. 

In the physical scheme one has to express SCETi form factor £ w through the physical form factor 
/+(0) = f M- ;| _ { 

Z*=CA5fo-CAD CB1 *£ 1 > ( 33 ) 

where C A0 and C B1 are the hard coefficient functions which appear in matching of the scalar heavy-light 
QCD current to the operators (|2l)l |2"T|) : 

q b = C A0 * J A0 + C B1 * J m + 0(l/m b ) (34) 

Inserting equation (|33|) into ([31]) we obtain 

oti = fo v t /C A0 * cpn + ^ * (U - Vl C B1 /C A0 ) * £ B1 , (35) 

Comparing this expression with equation ()10[) we find 

Vi(u) = v t (u)/C A0 , (36) 
Ti(u, z) = ti{u, z) - Vl (u) C B1 (z)/C A0 . (37) 

These expressions define precisely the coefficient functions Vi and T, in the physical scheme through the 
matching coefficient Vi and ti of the effective operators 0\ 2 . Introducing perturbative expansions for the 
coefficient functions: 

C AO = 1 + ^C AO LO + ..., (38) 

C+ 1 — Cfo + ~^~Cnlo + —j (39) 
Vi (u)=vt°{u) + ^v? L °{u) + ... , (40) 

i«(u, zj = t LO (u, z) + f^f LO ( M , z) + ... . (41) 
we obtained for the functions Vi and Tj in (fT2")l and (fT3")) 

(42) 

^°C^o , (43) 
tf°C&) , (44) 
-^°C B \ )-C B hvr°- (45) 

One can observe that subtraction terms —vf°C AQ LO in equation (|4*3")> . —Vi°C B Q and — wf° C B1 LO in 
(HU) and (|4"5)) can cancel the contributions of the " factorizable" diagrams which can be considered as 
appropriate product of the two matrix elements: 

<(7r7r)JO? i2 |B)^~<7r|?r g |0) (ir\^Tb\B) (46) 

If it is fulfilled, then such diagrams can be ignored in the calculations of the coefficient functions 
and T NL °. 



v LO 




-NLO 
i 




rpLO 


= (tf° 


iNLO 


= (t¥ LO - 
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Figure 1: Graphical representation of the QCD matrix elements which we use to compute Tj. The right 
graph corresponds to the QCD diagrams with emission of collinear gluon by collinear outgoing quark. 
Such graphs include contribution from operator j( A0 > with quark-gluon vertex from the leading order 
SCET lagrangian. 



In our paper we shall obtain the coefficient functions Ti computing matrix elements with quarks and 
gluons. For that purpose we define the perturbative analogs of the discussed form factor ^ 1 and LCDA 
<y3 w . To define SCETi form factor let us consider as external state a hard-collinear quark and a gluon 
with momenta pi and P2 respectively: 



pi = z p+ (pin)- +pi±, 
p 2 = z p + (p 2 n)- +p 2 j_, 

where p — rrib^, z = 1 — z and the other components are (n -pi, 2 ) ~ A, |pi,2_i_| 
for the hard-collinear momenta. Then we define 



(47) 
(48) 



\J Ami, as it necessary 



( q(pi)g(P2) I J {B1) (s)\ h v ) = m b f dz e is{2z-i){v.n)^ B e{z) _ 



(49) 



The factor £g denotes the relevant combination of the quark spinors and gluon polarization vector: 
Cb = -^x<l4g{\ — 75) ft»J3- We assume that the final gluon is transversely polarized with the polarization 
vector e g , the color indices are not shown for simplicity. Performing simple tree level calculation one finds 



3L0 



(r)=S(z-r). 



(50) 



In order to define the perturbative analog of the pion LCDA we consider quark-antiquark state with 
collinear momenta 

(51) 



Pi = W' + PlX, P2 = U P' + P2±, P' = m b 2 ■ 



Theifl 

(p' lP ' 2 \u(\n) i 75 d(-Xn)\ 0) = -1 f P (p'n) / dxe lX{2x ~ 1){p n) VP (x), (52) 

J 

where we introduced notation fp =i urjl^/su/mb- Again, from the leading order calculation one obtains 

ip^°(u)=5(u-x). (53) 

In order to calculate coefficient functions Ti we introduce the matrix elements describing the decay 
of the &-quark into three quarks and gluon 



(pi P ' 2 , Pl P 2 \ClOi + C 2 2 \ Mnf - {P1,P2 P1P2 \T {tl * J [M \ 4"c } I 6») nf = JpZb 



(54) 



and parameterized by form factors xf respectively. By the subscript "nf we indicate that we exclude 
the factorizable diagrams (|46[) which, as expected, cancel in the transition to the physical scheme. 



4 Symbols q and h v denote the quark spinors in this formula. 

5 For simplicity, we do not introduce here the collinear SCET fields following standard QCD notation. 
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a b 

Figure 2: Leading order diagrams for the coefficient functions Tj 



The subtraction term in the left side of eq. ([51|) represents the admixture of the operator J^ A0 ^ with 
one insertion of the interaction vertex £1™* from the LO hard-collinear SCET lagrangian [24J [25] . Such 
contribution describes the emission of collinear gluon from collinear quark and are present only in the 
diagrams with topology, given in Fig[T] Practically this subtraction can be easily done by the substitution 



ln¥g 2 (pn) 



(55) 



where u(pi) is the wave function of collinear quark in full QCD and u n denotes the hard-collinear spinor 
in SCET. 

From the factorization we expect that 



1 f 1 ~ 

du ipp(u) / dz ti(u, z) 0(z), 

Jo 



(56) 



where by tilde we denote coefficient functions of nonfactorizable diagrams. 

Consider, as example, calculation of '. Relevant tree diagrams at the leading order are shown in 
FigEl Two diagrams with emission of a gluon from the bottom lines represent the factorisable contribu- 
tions (|3S]) which cancel against —v[°C^q in ([J3|. 



Straightforward calculation gives (after Fierz transformation) 

Z? a = 0, 



i i ' r < I Ci±i 1 

D b = im b f P t B [- — - 



Comparing with ea. (|56p we obtain: 



du <p P °(u) / dz T 2 LU (u, z) 9 LO {z) = (i ml M B ) D b 



(57) 
(58) 

(59) 



Inserting the leading order expressions for the form factor 8 LO and LCDA ipp° we find the leading order 
hard coefficient functions: 



Ti(u,z) 



N r 1-U 



(60) 



As one can observe, LO results have no z— dependence. In the next section we use the same technique 
to compute the next-to-leading order corrections. 



3 Calculation of the coefficient functions Ti in the next-to- 
leading order 

Corresponding one-loop diagrams are shown in FiglHl Factorisable diagrams, in the sense of ((46|) . are 
not shown for simplicity. These are the diagrams where the external gluon is emitted from one of the 
bottom quark lines and the virtual gluon connects only the bottom (upper) quark lines but not upper and 
bottom. For the case of form factor ai, these diagrams naturally reproduce subtraction term Vi°C^\ 
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Dl-7 El-6 F5,6 G 1-5 Hl-6 

Figure 3: One loop diagrams which have to be computed in QCD. The crossed line denotes the emission 
of the outgoing gluon and the number gives the index of the corresponding diagram. The momentum 
flow is shown only for one line with the loop momentum I . We skip for simplicity the diagrams with the 
light quark wave function renormalization. 



in (|8gp and therefore cancel. But analogous situation for the a?2 is more involved because of the different 
Dirac structure of the operator vertex. The problem is that corresponding UV— divergent diagrams 
in dimensional regularization can not be represented exactly in the factorised form (J46J) because Fierz 
identities can not be used to regularized diagrams. Therefore one has to check the exact cancellation 
against v^C^lo a fter UV— renormalization. We have fond that in accordance with subtraction scheme, 
described below in the text, such cancellation is exact. Therefore we shall not discuss these diagrams 
further. 

From the factorization we expect that form factors x[ describing the matrix elements (|54[) of the 
renormalized QCD operators can be represented as a sum of three contributions: 

(*f ) NLO = rt°(x')*t» L °(x>,z>)*e L °(z') 

+ ^ LO (x') * tf°(x', z 1 ) * 9 LO (z') + <p P °(x') * tf°(x', z') * 9 NLO {z') (61) 

where fp LO and Q NLO denote the contributions of the renormalized matrix elements (|49|) and |52|) in 
the next-to-leading order. The three contributions in (|61[) can be associated with four integration regions 
in the loop integrals. The hard region ki ~ n%b, k 2 ~ mi provides contributions to the tf LO (x\z'), 
the collinear to p' must be associated with contributions to ipp LO (x'), the soft ki ~ A, k 2 ~ A 2 and 
collinear to p regions can be associated with the 9 NL (z 1 ). Substituting in (|61|) the explicit expressions 
for the ipp° , 6 LO (z') and tf°(x',z') we obtain 

. ( X J Uo + ^J ( '^fE^l + ^W^,, . (62) 

Inserting this expression into eq. (|45p and substituting C^q = —1 22, 26J we find for the NLO coefficient 
functions in physical scheme: 

Tr°( U , Z) = K ) NLQ + ^ jf 1 dx'^f^l + ^\[ dz>8^0 {z>) + yNLO. (63) 

Let us now discuss the calculation of different terms appearing in (|63[) . To perform the calculations 
of the diagrams in Fig[3] one has to introduce regularization for the ultraviolet (UV) and infrared (IR) 
divergencies. We shall use dimensional regularization with D = 4 — 2e to subtract UV— divergencies. 
To compute the UV— divergent subdiagrams of the four-fermion operators we use NDR— scheme with 
the anticommuting 75 matrix. Note that Fierz identities then can be used only for the renormalized 
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matrix elements in four dimensions. For the IR— divergencies we use regularization by off-shell external 
momenta. Such regularization makes possible to perform all manipulation with Dirac algebra for the 
UV— finite integrals in four dimensions. 

As an illustration, let us consider a contribution of some diagram Dx which can be represented in 
the following way 

D x = [ d D l uT x b v Jr 2 u, (64) 



where u, u,b v ,d are quark spinors of given flavor and matrices ri £g>r2 denote some momentum dependent 
expressions. Contraction of spinor indices can be organized in two different ways which correspond to 
the amplitudes a\ or a 2 - As an example below we consider the calculation of a 2 . The same technique 
also was used for a±. 

All graphs can be divided into two groups: UV— divergent and UV— finite. UV— divergent subgraphs, 
which appear in graphs G{1, 3, 5},H{1, 3, 6} represent usual divergencies of the QCD Green functions. 
They are removed by QCD Lagrangian counterterms. UV— subgraphs in diagrams A{4, 5, 6}, £?{4, 5, 6}, 
C{4, 5, 6}, -D{4, 5, 6} and ^{5, 6} describe renormalization of the four-fermion operators 0\^ 2 . Calculation 
of the corresponding UV— divergent integrals must to be performed in D = 4 — 2e. A typical expression 
for the integrand of UV— divergent graph can be written as 

^ b -^' N ^W\ +N 'W\ +N °W\ (65) 

where the ^-independent functions N^'" contain Dirac structures and spinors from the numerator, D[l] 
denotes the denominator, which behaves at large Euclidian momentum / as D[l] ~ (I 2 ) 3 - Such situation 
is usual for the UV— divergent graphs mentioned above, except only diagrams with quark self-energy 
subgraphs. Substituting (|65[) in ([64]) we obtain 



1 ^l l V 

where we introduced 



D x = N% v J[l»l v ] + JVf J[^] + A^o Jo, (66) 



[<l n l'-fyf ( 6? ) 

and similar for others integrals. Taking into account the behavior of the denominator at large momentum 
D[l] ~ (Z 2 ) 3 , it's clear that only J[l^l v ] is UV— divergent. The other three integrals can have only IR- 
divergencies, regulated by the off-shell momenta and therefore can be considered in D — 4. One can 
easily express the tensor integral J^plv] through scalar integrals: 



J[l^u] = g^Ji + + n tl n u )J 2 (68) 



with 



Ji = ^^AJ[l 2 ]-J[{l-n){l-n)]), (69) 
J 2 = -J—.((2-e)J[(l-n)(l-n)}- J[f}). (70) 

Both scalar integrals J[l 2 ] and J [(I ■ n)(l ■ n)] have UV— poles. But in the coefficient J 2 the poles cancel. 
Hence we must contract g^ v with Dirac structure N% in D = 4 — 2s and expand the obtained expression 
up to terms ~ e. The reduction of all one-loop Dirac structures to tree spinor combinations can be 
performed with the help of NDR prescription 

yyiy«(l _ 75 ) ® 7m7aii7a12 (1 - 75 ) = 4(4 - e )y (1 - 75 ) ® 7|t (l - 75 ) (71) 

Let us here make following observation. Computing the integrals J[l 2 ] and J[(l-n)(l-n)] we do not remove 
the IR— regularization because these integrals can also happen to be IR— divergent and we must avoid 
the mixing of UV— and IR— poles. However, we observed, that one can always choose some "convenient" 
momentum flow for which these integrals have only UV— divergences and free from /i?— singularities. In 
this case one can drop the IR-regularization making the calculations more simple. 
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Defining UV— pole of the coefficient J\ as 



pole part \J X \ = -Z uv (72) 



we write UV— divergent contribution for the Dx as 

D u x v = -Z uv lim N% u 9llv = im 2 b f P ^ B ^X col JX uv (73) 

where A co / is the color factor of the diagram and JXjjy represents some UV— divergent expression. In 
Appendix we provide explicit expressions for X co i and JXjjy for each diagram. 

The singular contributions D x v are removed by tree diagrams with operators renormalization con- 
stants and by QCD counterterms, see details in the Appendix. After that, performing Fierz transforma- 
tion, we obtaine trivial contributions from all the diagrams with topology ^45 — G5 and G{1, 3}. 

The similar calculation of the self-energy one-loop diagrams is much simple because they have only 
UV— poles and can be easily reduced to the tree Dirac structures. Let us mention, that we must consider 
also the diagrams with the one loop corrections to the wave functions which are not listed in the Fig|3j 
except for the heavy quark £&. We discuss these contributions in the Appendix. 

Finally, the finite expression for the diagram under consideration can be written as 

[°x]r = [K U 9^Ji] R + AT(rvn„ + n„n v ) J 2 + N?J[l^] + N J , (74) 

where [...]r denotes renormalized quantity. Remaining integrals have only IR— divergencies: collinear and 
soft. Calculation of such contributions is the same for all diagrams, with and without UV— subgraphs. 
Because we use off-shell regularization, we put D — 4 and perform projections on the pseudoscalar state 
with momentum p' : 

u a p u b a = H) /p^" W'-ftU , (75) 

and on the SCETi operator J B1 : 

d a a b c vp = d n t A ri_K [V^] 0a ®tt+ , (76) 

(77) 



[Vbi]/3q ~ 9 



13a 



where t are the standard color matrices satisfying Tr(t A t B ) = 

^5ab and dots stand for the irrelevant 
spin and color structures. It is convenient also to insert parametrisation for the vector integral 

J[l^ = \n^[{l-n)] + \n^J[{l-n)]. (78) 

After calculation of traces and contractions we arrive to the expression for [1?^]^ which can be written 
as a sum scalar integrals 

[D x ] R = im 2 b M B ^ [J x ] r , (79) 
[Jx\ R = [a 21 (e)J[l 2 ] + a 22 {e)J[{l ■ n)(l ■ n)]] R + b n J[{l ■ n)} + b 12 J[(l ■ n)} + c J . (80) 

with some coefficients a 2 i, bu and cq. The coefficients a 2 i, in front of UV— divergent integrals are computed 
in dimension D and therefore depend on e. 

The sum of the integrals [</x] fl gives the formula for the form factor xj : 

(*?)nlo = H (81) 
x 

Let us stress again, that expression with brackets [■■■]r in l|80p denotes the renormalized quantity. We 
simply have rewritten the coefficients J\^ 2 in terms of the corresponding integrals (|7T)|) . ([55)) and introduced 
e-dependent factors a 2 i which arise from the calculations in DR. Assume for simplicity that UV— divergent 
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integrals are free from IR— singularities. As we have discussed above, such situation can be realized for 
each diagram. Then contributions associated with these integrals are simply some finite expressions which 
do not depend on IR— regularization parameters. Our task now is to compute the remaining integrals 
J[{l-n)\, J[{l-n)] and J . 

Evaluation of these four-dimensional integrals can be performed with the technique known as expan- 
sion by regions, see for instance [27) . The dominant regions have been discussed above in the text. Hence, 
instead of one finite integral we obtain the sum of more simple but divergent integrals. According to gen- 
eral prescription, we use dimensional regularization in order to regularize the simplified integrals in each 
region. Therefore, in accordance with the dominant regions we have find following general decomposition 
oi[Jx] R 

\ J x] R = (Jx) hard + (Jx) 

coll—p + (Jx) soft ■ (82) 

Taking into account that IR— divergencies are related with collinear and soft regions we find: 

(Jx)hard = Me)J[Z 2 ] +a 22 (e)J[(l ■ n)(l ■ n)]] R + 

(buJKl ■ n)\ + b 12 J[(l ■ n)\ + co J ) hard , (83) 

(Jx) coll , so ft = (hiW-ny+hiW-ny + coJo)^,^. (84) 

The hard region contributions in DR have IR— poles instead of IR— logarithms as in the case of off-shell 
regularization. The collinear and soft contributions have only UV— divergencies (off-shell regularization 
works in IR— regions). But the sum of all terms must be finite because of cancellation between UV— 
and IR— poles. 

The contributions from the collinear and soft regions depend on the external off-shell momenta which 
we use as IR— regulators in the original integral Jx ■ Inserting decomposition (|82|) into (|8Tj) and then into 
([55]) we must recover the cancellation of the soft IR— scales. This is a good check of the factorization in 
the next-to-leading order. It is convenient to define the quantity: 

S = X! {( J x)coll-p + ( J x) coU - p > + (Jx) soft j 
X 

r- , -, / r 1 ,n NLO (T'\ 1 r 1 







dx ,Zr + dz'0 Nhu {z') \ . (85) 

1 - x' U . n J 



Then taking into account that ipp LO and Q NLO are defined as matrix elements of the renormalized 
operators ( the UV— poles are subtracted )0 we expect that the answer for S can be represented as a sum 
of UV— poles and scale independent constant: 

S = \z 2 + ~Z t + Z Q . (86) 

e e 

The poles, arising from the collinear and soft integrals in (|86p . must cancel against IR— poles appearing in 
the hard integrals (Jx)hard m (EU)- It is clear that the residues Z 12 can be related to the corresponding 
renormalization constant of the LCDA (pp and form factor 9. 

In order to obtain the finite terms Z one has to perform calculation of the collinear and soft integrals, 
matrix elements (|4^|) and (fB"2"]) and compute the sum. It is clear that both calculations overlap and this 
may be used to simplify derivation of the term Zq. The SCETi non-renormilized matrix elements must 
be computed in DR, as usually, with D-dimensional Dirac algebra. But the structure of numerators of the 
diagrams for the matrix elements are relatively simple and their reduction to the basic combinations can 
be evaluated without any special prescriptions. The calculation of the form factor (|4T)|) is more complicate 
in comparison with pure collinear pion LCDA ([52")) because certain diagrams generate UV— poles of second 
order. But corresponding diagrams always have very simple numerators. Performing reduction of Dirac 
algebra to the basic structures one can again rewrite expressions for the matrix elements in terms of the 
scalar integrals which are similar to those appearing from the collinear and soft regions in QCD diagrams. 



Note that one must apply the same IR— regularization for the matrix elements which define ipp LO and 8 NLO . 
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Now the coefficients in front of these integrals are e-dependent. We have found following representation 
for the bare NLO SCETi matrix elements: 

G±i ( f \ ,v n p lo (x') , i rl 



N c 



{(Jx) coll-p + (^x)coU-p' + (Jx)soft }+e^/x (87) 



X " X 



where Ix are some UV— divergent integrals. These integrals have only first order poles in e. All integrals 
with second order poles are absorbed into the first sum in rhs ([87)1 . Therefore taking into account that 
integrals Ix have coefficient e we obtain finite contribution from the second term in rhs ([57)1 . The factor 
e in the numerators appears, as a rule, from the reduction of D-dimensional structures in the diagrams 
to the basic factors £b and fp. Combining ([57)1 with ([85)) and taking into account UV— counteterterms 
for the SCETi operators we obtain representation ([86)) . It is clear that terms e Ix provide contributions 
to the constant Zq in ([55")) . 

Substitution ([2"T)) and ([55)) in formula for the coefficient function ([55)) gives: 

7f LO ( U , z) = £ (J x ) hard + ±Z 2 + \z x + Z Q + ^V(u). (88) 
x c 
This is our final working formula. It is convenient to rewrite the pole contributions as 

z, = (zo s/coll-p + {Zl)coll-p< (89) 

where (^i) s / co ;i_ p and (^i) co H-p' is associated with the form factor Q NLO (soft and collinear - 
p regions) and LCDA (pp LO (collinear-p' region) respectively. Calculation of the relevant integrals 
(details are considered in the Appendix) gives 

(Zdcau-p, =C F ~ (2 + In u), (90) 

V Nr. I U 



(91) 



1 1/ x / C+, 1 \ / 1 Cf , r ,/ , - , 1 r lnz C A z + lnzl\ 

^ + 7^)./^= (-if.) (-^f (l + 2 £ ln[ M /W) + -^--f — ]j 

We checked that above expressions are consistent with the evolution kernels of the pion LCDA [28j [29) 
and SCETi operator J B1 [T3l [24]. But note, that we haven't included the diagrams with the wave 
function renormalizations hence the poles ([90)) and ()9T)) are not exactly convolutions of the corresponding 
evolution kernels with the leading order coefficient functions. 

Calculation of the hard integrals ( Jx)i ia rd can ^ e done with standard technique and the results for 
each diagram are listed in the Appendix. The arising IR— poles cancel in the sum with ([90)) and |9T)) as 
it is expected from the factorization theorem. Resulting expressions for the coefficient functions can be 
written as 

a 2 : T 2 NLO (u, z) = ^f D (u, z) + ^-T ND (u, z), (92) 

ax : if LO (u, z) = ^-Td(u, z) + ^T ND (u, z), (93) 

where the subscripts D and ND can be understood as diagonal and non-diagonal contributions. Assum- 
ing for simplicity fj,R = /j,f = pb we obtaine: 

1 l/4-u 2u 2 u \ 
— lmJ c h, z) = — tt — 

it m ' 2 \ u (1-u-zf l-u-z) 

(l zu 2 \ 
+ k mu 

\u a-u- z? ) 



■\\\u + u\-- H 3 lnz , (94) 

(1 — u — z) I 



u u — z 
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^ m i \ 3 1 / o / o \ 8 1 / u 2 + 3« 2u 2 u u \ , 

Re T D u, z = - In (/i 2 /mg) + - + — + 3 + , ^ Inu 

u v ' u 2 1 u 2 (1-u-z) 1-u-z u 2 z/ 

(1 u u \ _ _ 

- 9 _ : 3 - - 9 r (1 - uz) ln(l -uz) 
2u 2 z u[l-u-zf 2u 2 [l-u-z)J y 

, / 3 1\ (2-3*) , 

■ ln(l — uz) — — H Inn - _ M nz 



uuzz \2u uz J 2uz 

z(z — 3u) , ln 2 u ln 2 u lnz,, 
— ^-2 lnz - -^r- + -7^- + -p- (Inu - lnu) 



2(l-u-z) 2 2« 2u u 

Li 2 (u) - Li 2 (u) z u 2 uz T ,_ . 

+ — - ra /(«, *) + -=-/(«, ^) , (95) 

u (1 — u — z) u 



1 C A / u l + 2u 

■ Im T/v D (u,z) = 



7T ' ' ' 2 V 1 — u — z 2u 



In u + tt — 1 hi 

l(l-u-z) 2 J 



u — z 



z , / 1 u 2 u \ , 

+ — lnz+ o + lnz 

uz \^u [1-u-z) u-zj / 

jl/5 — u 2u 2 u ^ 

+ C F - 



2 ^ u (1-u-z) 1-u-z J 

1 u 2 z \ wz,_z 

— o In u + — In u — In z 

u [1-u-z) ) u[u- z) uz 



[i-u- zy 



u u — z 



lnz (96) 



Re T ND [u, z) = C F T CF [u, z)+^ T C a(u, z), (97) 

T C f(u, z) = -L In 2 (u 2 /m 2 ) + 1 ^ + Inn - ^ In ( M 2 /m 2 ) 

27 7T 2 (4 - 3 u - 5 u z - 4 u 2 z) ( 1 2 \ 
H — r H = - ( — H = lnu 

2 u 24 u u z \2u uz J 

( 2 + u u 2 u 1 \ , 

+ — + 7 72 + TTTi 7 + — lnu 

\ u (1-u-z) 2 2 (1-u-z) uz J 

3 (3- 2z-uz) , 2 (1 - uz) , , 
tttt^ -In z H _ _ ln(l - u z) 

2uz uuzz 

1 / 3 2u 2 u 1 3-5u\ , 

+ - - + 7 72 + 1 2 + = — ln * 

2 \ u (1 — u — z) 1 — u — z uz^ uuz / 

+ / 1 l (2-u)u 2 u(u 2 + 2u) 



2uuz 2 u(l-u-z) 2 2u 2 (1-u-z) 

6u - 3 - 4u 2 1 \ , . 

+ 7T^2 — In 1-uz 

2uu z z uz J 
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Inn / _ I + uz 
— — — (1 + u) m u H — In z — In z 

u \ z 



ln 2 u ln 2 u (2 - z) 



u 2 u 2uz 

z — 2 u z 



In 2 z- 



(1 + u) 



In it In , 



M Z 



-Li 2 (u) 



1 y . /— \ 1 - u (3 - z) 

Li 2 (u) Li 2 (z) 



iia z 



(1+Mz) 

— — Li 2 (it z) 



— ^ - 1 ] Li 2 (Sz) + \hi 2 {z) 
uuz ) u 



u 2 z uz 

-$Hu, z) + —I(u, z) , 

(1 — u — z) u 



(98) 



T C a(u,z) 



In z 



In (/i 2 /m 2 ) 



9 7T 2 (1-it (3-z)) 
w 6uuz 



lnu 3 (2-^) 



zt z 

UUZZ 



2iiz 



In z 



2 m 1 — u — z 



In it 



, , , (1 — uz) (z — u (1 — z z)) , , 
ln(l -uz) + - z4p — rr^- ln(l - u z) 



uu (1 — u — z) z z 
u 1 — 2 M 



1 

2~S 

- In it 



2u I — u — z 
1 



lnz 



ln 2 z- 



ln 2 u In 2 it 



2u 



-ln 2 z 



■^■^ Inz + mz H — (in z — u m z) 

uz I u 



I 1 +UZ IT', 

1 — I Li 2 (u) 



1 



uuz 



Li 2 (u) 



1 T . , . (1 — 3tt + « z) T . . . (l + wz) T . , . 
+ -Li 2 (z) + i — i Li 2 (z) + v _ _ Li 2 (u z) 

u uuz uz 



(99) 



Li 2 (uz) - -Li 2 (z) + I{u,z) — u I{u,z) 

uuz ) u 1 — u — z 



(100) 



where we introduced convenient real function I(u, z) which reads 



I(u,z) = 



ln(w/z 

u — z 
1 



ln(uz) — ln(z u) — ln(u + z) 



ln(w/z) 



In 



u + z 



2Li 2 (z) + Li 2 



z 



u + z 



Li 2 



z 



(z <-> It) 



(101) 



Our results are in agreement with the kernels which have been earlier presented in the paper [15j but 
computed using different technical approach ( dimensional regularization with evanescent operators 
Let us also observe that the function I{u,z) introduced in eq. (|10ip naturally appears in calculations 
of the diagrams involving massive propagator of heavy quark. The similar structure have been also 



7 We are grateful to S. Jager and M. Beneke for the correspondence which helped us to fix a mistake in the expression 
for one diagram 
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introduced in [15] and denoted as F(z, u). Let us stress once more, that I(u, z) is a real function as one 
can easily see from its definition IjlOljl . We have found that 



I(u,z) = Be[F(l-z,l-u)]. (102) 

z — u 

We did not find any transformation to prove this equivalence analytically and checked it numerically for 
the several arbitrary values of the arguments. 

Analytical expressions ([94| - (|100p for the coefficient functions T/^; represent the main technical 
results of our paper. 



4 Numerical estimates of B —> tttt branching fractions 

In this section we perform the numerical analysis of the branching fractions including next-to-leading 
corrections to the hard and jet coefficient functions. The main contribution to the branchings originate, 
obviously, from the real parts of amplitudes a% t 2- We neglect in our estimates by electroweak penguins 
&i,EW but include QCD penguins a± ,c , sec cq.(Ql, in the form presented in |9J. 

Consider first some important details in calculation at the NLO approximation. General formula 
reads: 

oti = fa\ du Vi(u, n R , [i h )ip n (u, [i h )+ (103) 
Jo 



du <p„(u, hh) I dz Ti(u, z, (i R , fih) £tt 0*,Mfc) , (104) 
Jo 

/>oo /> 1 

£%(z,fMh) = /bItt duo I dx (j> B {u,liF)J(z,x,u,iih,ii F ) (p^(x,fi F ). (105) 



where we have shown explicitly the scale dependence. To estimate the values of these amplitudes we use 
the following numerical input. For the coefficient functions Ci=i,2 (/■<«) in the effective Hamiltonian ([3]) 
we employ the NLO results at [i R = rrib obtained in [T7] (A^£>i?-scheme, NLO approximation): 

C? LO (m b ) = 1.075, Cf i0 (m 6 ) = -0.170. (106) 

where m& denotes 6-meson pole mass rrib = 4.8 GeV. For two others scales we accept the following values: 

Hh = rn b , fip = fjt hc = 1.5 GeV. (107) 

For simplicity, the uncertainties in the scales setting will be ignored in our analysis. Then we obtained 

<*i=fo / duVi(u,mb)<Pir(u i rn,b)+ / du tp^ (u, m&) / dz T t (u, z,m b ) £f (z,m b ) , (108) 
Jo Jo Jo 

To compute the form factor ^ 1 (z,mf ) ) we must perform evolution from scale m b to the scale /i^ c : 

^(z,m b ) = / dz'U B (z,z' 7 m b ,fi hc ) ^(z',fi hc ), (109) 
Jo 

where the evolution operator Ub is derived from the solution of the evolution equation [TBI |2~4"] : 

-f- g(z, /i) = - / dz [V(z, z') - 6(z - 2')r c „ sp (a s ) In -^}tf(z', //), (HO) 
dhi^i J m b 

The explicit expression for the evolution kernels in our notation are the same as in [24], see Eq's (46,47). 
Let us write the solution for Ub in the form 

U B = U L B L { Zl z') + U^ LL (z, z') + (Ill) 
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where Ug L and Ug LL can be understood as "leading log" and "next-to-leading log" approximations. 
Both terms are needed to perform complete calculation at next-to-leading order. The evolution of the 
form factor is combination of the two effects: summation of the so-called Sudakov double logarithms 
associated with cusp-anomalous dimension r cMsp and usual light-cone logarithms associated with non- 
local evolution kernel V(z, z') (or non-local anomalous dimension). At present T cusp (as) and V(z, z') are 
known at two- and one-loop accuracy respectively. This is enough for summation of leading logarithms 
in Ug L but in order to compute Ug LL one has to know r cMsp at three loop and V(z, z') at two loop 
accuracy . Because these quantities are uncalculated, we can't perform corresponding evolution. So we 
just neglect by this effect in our calculation. 

Form factor (z, fi^c) can be computed systematically order by order using factorization approach 
. Corresponding jet function has been computed at the next-to-leading accuracy in Q21 EH 03] . Fixing 
the factorization scale at this step to be equal to jip = jih c we escape the summation of the large 
logarithms in jet function and therefore we have to provide the LCDA's (p v (x, fihc) and 4>b(w, Hhc) at this 
normalization point. Assuming decomposition — (£^) LO + ) NLO we obtaine: 



MSA 



du >»,,) / dz T t {u,z)U B {z,z') * [(^) LO (z') + ^-(^) NLO (z') 



= / du^(u)T t LU (u) / dz U^(z,z>)*(tf) LV (z>) + 



Zir 



du <p n (u) Tt°{u) I dz U£ L (z, z') * 



du <Pir(u) 



dz 



U^(z,z')*(^) LO (z') + 



(112) 



(113) 



where dots denote the neglected logarithms associated with Ug LL and higher order terms 0(a 2 s ). For 
simplicity by asterisk we denoted the integration with respect to z' and skip obvious scale dependence. 
In the first two lines of (|113p we used that T? 1 does not depend on fraction z. Then one can perform 
integration over external variable z 



f 

Jo 



dz U^ L (z,z') = U LL (z') 



(114) 



that simplifies the evolution. In the numerical calculations we have used for U LL (z') simple approximation 
that was found in [13]. Therefore we obtained 



dz U^ L (z,z')*(^) BLO (z') 
where we used standard notation 



N c K F X B m b 



-<t>B (u,Hhc 



— ^(z,fi hc )U LL (z), (115) 



(116) 



Assuming the following ansatz for the pion DA amplitude 

Pir(u, fihc) = 6w(l -u) + al(^ hc )Cl /2 {2u - 1) 
one obtains (useful technical details can be fond in [13] ): 



aUVhc)C 3 4 /2 (2u-l), 



/ — Mz)U LL {z) = 2.72 (1 + al 
Jo z 



al). 



(117) 



(118) 



The similar calculation for the next-to-leading contribution gives 



dz Ub (z,z) * 



CO 



d Z ' u LL (z')^r L °{ Z '^ hc 



(119) 
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otsfBU/Kp I dz'U LL {z') I duj I dx 4>b{uj) -^Jnlo(z' ,x,uj) tp n (x) (120) 



2tt 



-nas 



N c K F X B m b 
[3.0+ 10. 10aJ + 16.60 a J] ), 

where we introduced 

(L) = Xb / </>s(a;) In 



J ( (! + 4 + a I) {L 2 ) ~ I 3 - 93 + 8 - 15 «2 + 10.05 aj] (L) 



m b ui 

Vic 



m b uj 

tic 



(121) 
(122) 

(123) 



The last term in eq. (|113p with the convolution of the next-to-leading coefficient function can be repre- 
sented in following form: 



-irasinhc 



N c Kp Xb mb J 27r 



m,n=0,2,4 



(124) 
(125) 



where the moments £, 



a 

i 



1 ±mn 
I n 



a 



»±i 



A/'. 



f 1 



D(ND) 



duQuuCH 2 {u-u) I dz T% ( L ° D) (u, z) dz'U% L (z, z')6z C^ /2 (z - z) 



(126) 
(127) 



have been computed numerically. If the evolution is switched off our results are in agreement with the 
moments computed in [15], see equations (54,55). Performing the evolution and with the scale fixed as 
described above we obtained: 



= 73.64- 



tj° = 27.94- 

+40 
b n 



-12.74 



13.85i7r 
21.22i7r 
f 22.61Z7T 



j-02 
D 



t 



22 
D 
f 42 . 



66.80- 
23.86- 
-9.11 



16.70i7r 
28.46wr 
-31.37wr 



= 66.85 



<24 = 24.44 

j.44 



17.74m 
32. 10m 
7.71 + 34.67i7r 



+ oo _ 

ND — 
±20 _ 
l ND — 
+40 _ 
ND — 



-7.34- 

201.69 

371.10 



14.24i?r 

- 39.87i7r 

- 51.0H7T 



L ND 



= 33.15- 17.43i7r 



t 2 ND = -180.99- 39.19i7r 

f 42 
ND 



t°^ D = 69.63 



-347.00 - 49.66i7r 



.24 
b ND 
+44 
ND 



18.21WT 

= -164.74- 39.28i?r 
= -335.62 - 49.63i7r 



(128) 



(129) 



In order to compute branching fractions we used the set of input parameters given in the table below. 
The values for the coefficient functions 63,4,6, s e // are taken also from [T7] . 

Table 1. Numerical values of the phenomenological parameters used to compute branching fractions. 



u 


131 MeV 


fm 


200±30 MeV 


a£(1.5GeV) 


0.25±0.15 


A B (1.5GeV) 


0.35±.15 GeV 


<(1.5GeV) 


-0.05±0.15 


Jo 


0.28±0.05 



\V ub \x 10 3 


\V cb \x 10 3 


7 


TB, PS 


C 3 (m b ) 


C'i(m b ) 


C 6 (m b ) 


C 8 eff(m b ) 




41.4 


62° 


1.53 


0.012 


-0.030 


-0.035 


-0.143 



m B 


5.28GeV 


pole 
m b 


4.8GeV 


m pole 


1.8GeV 



For the pion LCDA we use a simple model with two Gegenbauer moments a\ 4 (|117[) . Our estimates of 
the moments based on the results obtained in [3TJ [35] ■ The evolution of these moments from initial scale 
fj,hc to scale m b have been computed with the next-to-leading logarithmic accuracy for the leading order 
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contribution (the first term in (|113[> ). To perform this two- loop evolution we have used the analytical 
results derived in [33] , 

For the B-meson LCDA we accept a simple model with exponential behavior which is very popular 
in phcnomenological applications 

<Ab(w) = -T2- exp (-u/X B ) ■ 
Then one can easily calculate the moments introduced in eq. (|123p : 



(L) = In 



rribXi 
Ac 



7B , <L 2 )=ln 2 



TOfcAi 



IE 



TV" 



where je = 0.577. 



With the given above central value of Xb one obtains 
(L) = -0.87, (L 2 ) = 2.4 



For QCD running coupling as we use the two loop approximation with QCD scale 1^q CD — 225MeV. 
Recall, that in our numerical estimates of the branching fractions we neglect EW-penguins contributions 
but include QCD penguins in the NLO approximation as given in [5]. As it was observed in those papers, 
the values of the pion branching fractions are very sensitive to the product | V ub \ fo ■ Corresponding 
value can be estimated from semileptonic decay B — > ixlv assuming monotonic behavior of the form 
factor f + (q 2 ) : 



dT(B° -► nlv) 
df 



Gp 
24tt 2 



\V ub \ 2 \f + (q 2 )\' pi 



Gp 



24tt 2 



\v ub nu( q 2 )\ z P i> 



Gp 
24tt 2 



\V ub \ 2 \f Q (0)\ pi 



Using results obtained by BABAR collaboration in 

2 dT(B° -> nlv) 



TB 



dq A 







dq 2 



for the lowest bin in q 2 < 8GeV 2 : 
ABr(B° -> nlv) = 0.21 ± 0.13 



one obtains 



10 3 | Kb | jo < 



1 ABr(B° -> nlv) 
W 1 Jo 8 dc ?Pl 



x 10 J 



n 79 +0.20 
u -' z -0.27 



1.0), 



where in brackets we show the product of the central values \V u b\ and fo from the Table 1. In order to 
satisfy this requirement (at least for upper bound) we accept following values for \V u b\ and fo : 



\V ub \ = 0.0038, f = 0.23, with 10 3 \V ub \ f - 0.87 



(130) 



As one can see from Table 1, such choice of the \ V ub \ and fo corresponds to the lowest possible value of fo 
within indicated uncertainties. First, we compute two largest branching fractions as a functions of four 
parameters Xb, anc ^ a \ ■ These results for the tree level dominant branchings Br(B — > tt~tt + ) and 
Br(B — > 7T~7r°) and corresponding values for 0:1,2 are presented in FigHJ We show all solutions which 
describe the experimental points changing the parameters inside the intervals indicated in the Table 1. 
As one can observe, there exist many possible solutions that demonstrate large ambiguity due to badly 
known mesons parameters. For instance, we reproduce the experimental values 



10 6 Br(B 
10 6 Br(B 



7T-7T+) = 5.1 (exp: 5.0 ± 0.4), 
tt-jt ) = 5.51 (exp: 5.5 ±0.6), 



(131) 
(132) 



with fs = 0.23, Xb = 0.23, a\ = 0.3 and a\ — —0.07. Corresponding amplitudes 01,2 have following 
numerical structure at this point: 



oi// = [1.04 + 0.012i] v + (-0.030)^0^0 
= 0.96-0.019i, 



a-ilh = [0.035 - 0.077i] y 
= 0.38- 0.020i, 



-0.020)^0^0 + (-0.035 - 0.03U) tnlo ,jlo 
(0.19) tlo ,jlo + ( 0.13) TiO , JNLO + ( 0.028 + 0.060i) T 



-■JV LO *jLO 



(133) 
(134) 
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Figure 4: Branching fractions and and amplitudes 0:1,2 = ai^/Zo as functions of hadron input parameters 
/b(0.02), Ab(0.06), aj(0.06) and aj(0.06). The numbers in the brackets give the value of step in the 
numerical calculations. The large points in the upper plots correspond to the choice fg = 0.23, = 
0.23, = 0.30 and a\ = —0.07. The values of the on and branching fractions which lie outside of 
experimental interval are not shown. 



where for convenience we presented the answers normalized to l//o- Let us briefly comment these results. 
The real part of the ol\ is clearly dominated by the vertex term, the corrections from the hard spectator 
scattering are about 5% in absolute size. As one can see from (|133[) the radiative corrections ( indicated 
as TL o t jNLo and t nlo^,jlo) numerically quite large with respect to LO term t lo*jlo. The relatively 
large value of T NLO * J LO contribution is explained by large value of the Wilson coefficient C% with 
respect to C2. 

For the amplitude a<i the situation for the real part is different. The vertex contribution is very 
small due to the compensation between tree and NLO contributions [9]. Therefore the dominant term 
arises from the hard spectator scattering part of the amplitude. The structure of the NLO terms here 
is also different: bulk of the radiative correction is due to NLO jet function. The NLO hard spectator 
scattering is approximately four times smaller. Hence obvious conclusion is that the total (jet+hard) 
NLO contribution is very important for Re[o2] and almost negligible for Re[ai] respectively. 

The important result of the our calculation is the imaginary part of both amplitudes 0:1,2 • Its value, 
in comparison with the value of imaginary part from the vertex contribution V, is quite large and has 
opposite sign. Therefore the resulting imaginary parts are significantly modified. Such changes may 
produce sizable effect for the CP-violating asymmetries and therefore have to be taken into account in 
phenomenological analysis. 

The smallness of the 02 provides small value of the third branching 10 6 Br(B — > 7t 7t°). Its value 
always remains considerably smaller than the experimental value, see two bottom plots in Fig ((4|. With 
the amplitude a-i from eq. (|134p we obtained 

10 6 Br(B -> 7T°7T°) = 0.45 (exp: 1.45 ± 0.29) 

that is three times smaller then the experimental value. Of course, there exist ambiguities not only due 
to hadronic input parameters but also in the scale setting, in quark masses and weak parameters. Such 
uncertainties have been already estimated in [15] and they are quite large. 

On the other side it's possible to suppose that realistic explanation of the small theoretical value of 
the 7T°7T° branching can be explained by relatively large contributions of the power corrections which 
have been ignored in present calculations. The key observation is that B —> 7t 7t amplitude has small 
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absolute value (~ as) due to cancellation among the vertex contributions. Then it makes possible that 
preasymptotic effects from the power corrections are very important especially for this case. In papers 
[9]fl5] the model for some contributions of the power corrections have been already introduced to estimate 
their effect. In particular, the so-called "chiral enhanced" twist-3 contributions have been considered as 
a source of dominant effect. In |15j such correction strongly enhances the absolute value of a 2 compare 
to perturbative contribution. 

Similar situation, with small leading power term and large power suppressed corrections occurs in 
hard exclusive processes. Such scenario, as expected, is realized for the pion form factor in large Q 2 limit. 
In that situation leading twist perturbative contribution is small ~ 0(as) and power correction may 
even dominate in quite large accessible range of Q 2 . For the detailed discussion of this question we refer 
to [34j where light-cone sum rule approach have been used for analysis of the power behavior in Q 2 . The 
other interesting observation, which was made in 34], is related to the "chiral enhanced" contributions. 
For the pion form factor such corrections naively could provide very strong effect as it was observed first 
in |35j . But in sum rule calculations it was found that such contributions turn out to be small. It might 
be understood, that the value of such corrections is overestimated if one uses the simple model with a 
cutoff of the momentum fraction to ensure convergence of the convolution integrals. If this true, then for 
description of B — > mr decays one needs, probably, a different model of the power corrections than one 
used up to now. The detailed discussion of this question lies beyond the scope of present consideration 
and we refer to recent works dedicated to this subject [36] . 



5 Conclusions 

We have presented the independent calculation of the next-to-leading order corrections to the graphical 
tree amplitudes in B — > ttt: decays. Our analytical expressions are in agreement with results obtained 
in [15] using a different technical approach. The obtained results have been used for the numerical 
estimates of the branching fractions in BBNS approach [9]. We have found that total ( hard plus jet 
) next-to-leading correction is relatively small for the real part of the a± decay amplitude and provide 
large contribution to the real part of a 2 . In the last case the dominant effect originate from the next-to- 
leading contribution of the jet function. But the imaginary part which is generated by the hard coefficient 
function of the hard spectator scattering term is quite large and therefore can provide sizable contribution 
to the CP-violating asymmetries. Our estimates of the branching fractions allows to make conclusion 
about existence of sizable effect from power suppressed contributions, especially for branching B — > 7r°7r° 
which dominated small a 2 amplitude. 
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Appendix 

1. Structure of the divergent contributions of the Feynmann integrals 

Here we briefly discuss the structure of different divergent contributions and provide results for the 
singular parts of the UV— and IR— integrals. 

The UV— divergencies arising in the NLO diagrams are removed by the counterterms for the four-quark 
operators, renormalization constants of the wave functions and QCD counterterms. Renormalization of 
the four fermion operators 0\_ 2 is given by: 

O x = Z n Z^ 2 O b r e + Z 12 Z^ 2 O h 2 are , (135) 
2 = Z 12 Z7 2 O b r e + Z 22 Z7 2 O h 2 are (136) 
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where 2x2 leading order matrix Z reads (see for instance |17| ) 

2C F - 6(C F - C A /2) 



„ 1 a, 
Z= 1 

e 47r 



-3 

2C F - 6(C F - C A /2) 



(137) 



and Z$ is the renormalization constant of the quark field in the MS-scheme. 

v e 47r 



(138) 



In addition, matrix elements have to be multiplied on the corresponding renormalization constant of the 
wave functions external particles. For the quark wave function such renormalization factor is defined by 



1+ t- 



dp 



(139) 



p=p 



where we assume off-shell IR— regularization and denotes one-particle irreducible self-energy graphs. 
Such definition introduces, except UV— pole, the finite term: 



In 4 



Zfin 



(140) 



The pole part cancel UV— divergencies arising in diagrams for the matrix clement. The finite term must 
be included in the matching. But the same definition (|139p is used for the renormalization of the quark 
wave functions in the effective theory. One obtains similar contribution: 



1 



(Z g - 1) 



47T 



-c t 



hi 



J fin 



(141) 



which is different from (| 140[) only by renormalization scale /.ip. Hence, the difference of two expressions 
(|140p and (| 141 [) which defines the contribution to the hard coefficient function from such terms is pro- 
portional to [i F ). We put /in = [i F that allows to avoid consideration of such terms. The same 
arguments can be repeated for the gluon wave functions. But situation with heavy 6-quark is different. 
The HQET wave function of the effective filed h v is renormalized by factor Zh 



1 + i 



d(vk) 



(vk)=0 



where k denotes residual momentum. Then one obtains: 



2 In 

e 



1 



as 



-{Z h -l) = -^C F [ — -2\n 



■in 



2(vk) 

2(vk) 
jj, F 



31n[^ 



Hence the corresponding contribution to the hard coefficient function ( for /i^ = /i F ) 



as 
2tt 



C F 



i±l 



N c u 



2 nib 



(142) 

(143) 
(144) 

(145) 



After these remarks let us provide list of the singular contributions for the integrals (Jx)uv defined 
in (|73[) and for the soft and collinear integrals (|8"2")l . ([84]) which appear in graphs presented in Fig[3J For 
simplicity, we shall indicate below index IR for the soft and collinear integrals and define (X — Al, A2, ...) 



pole terms [(J X ) hard \ = -pole terms {Jx) soft . coa 



— X rn i J X 



(146) 



where X co \ denotes color factor as in (j73|) . The explicit expressions for X co \ are listed in Table 1. 
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Tabid. The color factors of different graphs in Fig. [3] 



Color factors X co i 


diagram index X 


Qi , c i±t i r c A \ 


Al, 56, CI, L>6 


Ci , Cj±i (~i 
2N C ' N c ° F 


A6, C6, £>1 


C i± i /y-, CU\ 


A2, C2, Z32, E6, F6, HI 




A4, A7, 54, C7, £>4, C4, M, E2, HQ 


C i± i C A 
N c 2 


A3, C3, L>3, 53, 53 



As usually we use notation u = 1 — u and MS-scheme for subtractions. For simplicity we put fx = nib- 
Then 



JA6 



uv 



JC6 



uv 



JE6 



uv 



JF6 



uv 



J HI 



uv 



-J HQ 



uv 



.-L JH3 uv --^- £ , 



JAA UV = JC4 UV =~, JB4 UV = JDA UV = JB6 UV = JD6 UV = 4, 

2e £ eu 



JA1 



IR 



U \ £ 



' '' 2 + - (1 + m - \n\zu}) ) , JA2 IR 



1 1 

e u 



1 z In z 

2 F~ 



1 1 



1 1 



JA3/ fi = -- (1 + \nz) , Jv44 /i? = — + - (2 + in - ]xi[zu}) , JA7 



7i? 



1 1 

£ 2S' 



J54 /7V = + - (-2 - ivr + ln[zu\) , J56 /i? = 4 ( -4 - - (2 + in - \n[zu}) 
e e a 



1 



1 



JC2/R = — — In 



JC3m — 



2s 1 

1 

2 



JC4r H = - 1 - 



J5>1 



1 



IR 



1 



1 — z« 



1 



In 



in — In u 
1 — zu 



1 — zu 

uz 

Z + 2uz 
uz 



■ 111 



z 



In z + 
1 



£~ £ 
1 — Z?i 

1 - z 
1 — zu 



zu 

1 



In [1 



, JC6 IR = - — + -(l-lnu; 

U \ 2£^ 



2£ 2 



1 / 1 



£ 2 V2ft 



-1 - 



1 



1 



- -— 7 + -lnS , JD2 IR = - — + - lnu 



t7r In it 

2-0 



2£ 2 

(1 - 2uz) 



lnz + 



1 

£ 

1 — UZ 

zu 

(1 — uz) 
uz 



In 



1 — MZ 



ln[l — uz] 



JD4 TR = - -1 



J56 Jfl 



1 — uz 1 — uz 
m 



1 



1 1 

£ U 
1 1 



UZ 

1 

-- + - 

£ z £ 

lnu 

u 



„ 1 In u 
JE2 = , 

£ U 



1 / 1 



In? 



+ In [it z] 



1 



JHl IR — --, JH3 IR = — I — + - - + In [zw] 



£ It 



u \2£ 
1 /3 



' +i(l-lnz) 



(147) 
(148) 

(149) 
(150) 
(151) 

(152) 
(153) 
(154) 
(155) 
(156) 
(157) 
(158) 
(159) 
(160) 



Performing analysis of the main regions which contribute to the leading power accuracy we find that 
many IR— contributions cancel in certain combinations of diagrams. Important that such cancellation 
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can be observed before computing of the integrals and therefore provides a good check for the interme- 
diate calculations. Taking into account the color factors of diagrams we have found following IR— finite 
combinations: 



(JA1 + JB6 + JC6 + JD1) IR = 0, (JA4 + JB4 + JC2 + JCA + JD2 - 
(JA1 + JB6 - JE3 + JE6 + J HI - JH3) IR = 0. 



JDA) IR = 0, (161) 
(162) 



The remaining IR— contributions can be associated with the matrix elements of the SCET operators 
which define LCDA (fp LO and from factor 8 NLO . To obtain results for the ^0,1,2 introduced in (|137| we 
have computed following expressions: 



e 



a 



i±i 



coll—p' 



C«±i 



C F 
C f 



(JE2 + JE6 + JHl) coU _ p ,+ dx''??—)- 

Jo 1 x 



U(2 

e u 



In i 



u 



(163) 



s/coll—p 



C F (JA2 + JA7 + JF6) 



s/coll—p 



Ca 



(J A3 - JA2 - JC2 + JC3 - JD2 + JD3 - JF6) 



±1 



C F - 



1 



1 



u \ 2e 2 



Xhxz 

e z 



In z 

z 



CaI 
2 u 



s/coll—p 

1 z + In z 

e z 



dz'O 



inNLO 



1 



In z 

z 



In z 

z 



(164) 



The pole contribution in formula (|163|) can be interpreted as convolution of the evolution kernel with 
the leading order coefficient function T^: 



(ZJcou.p, + = V *tf° * 9 LO 



1 ( 3 , 
— — h in u 

u\2 



(165) 



where in the left side of eq. (|165p we introduced contribution from the quark field renormalization, denoted 
as Zjp. We added this term because the set diagrams in (|163[) doesn't have such contribution and therefore 
corresponding poles define non-trivial but not complete part of the evolution kernel V: 



V{ Xl u) = C F 



(x < u) 1 + 



1-u 









- 1 )1 




x — u J 



(166) 



where plus-prescription denotes: [F(x, u)] + = F(x, u) — S(x — u) J Q dx'F(x', u). The same consideration 
can be carried out for the poles in (|164[) . Because leading order coefficient functions (|6H)) are 

independent on the momentum fraction z', the singular (pole) part of the form factor 9 N (z') appears 
as integral J dz' 9 NLO (z') and can be understood as counterterm of the local operator J B1 (s = 0) in the 
effective theory. We have checked by independent calculation that expressions for the Z\p in (|164[) is in 
agreement with the renormalization of the local SCETi operator J B1 (s = 0) l|27p: 



(Zi)s/coii-p + Z^f, + -2^4 



= <p? J *T i 
C t±1 1 

N B U 



UV-pole — 

-i^(l + 2 £ ln[ M /(pn)]) 

£ 2 

z — 4 In z 



C F - 



4(1 -z) 



Ca In z 
~2~1 - z 



(167) 



where we again introduced contribution from the renormalization factors for coupling Z g and fields Z^^a 
which necessary for complete definition of the evolution kernel. Equation (|167[) is in agreement with 
known results for the evolution J 31 obtained in [T3l [24l. 
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2. Finite contributions of the hard integrals (J 



x ) hard 



In the second part of Appendix we present the finite contributions of the hard integrals (Jx)hard intro- 
duced in (El) and (El). We shall write 



finite terms [{Jx) hard ] = X col JX (168) 

As usually, we assume /j, = rrib in order to simplify the formulas. 

2 7ir 2 lnu (l + z)lnz (3z-2)lnz ln 2 (z u) in (1 - \n\u z}) 
u lZu u Zuz Zu z Zu u 

1 lnu / 3z , \ lnz z\n 2 z ( 1 zlnz\ , _ . 

JA2 ^—u + ^ + {-^ + zXnU )m + ^¥ - i7r {M + ^z-)> (170) 

Lb £j Lb \ Zj / Hi £> £j LI Aj \ Zj (J, Lb A/ I 

JAO _ 2 ln(l~") . Ml-z) ln(l - u) ln(l - z) ln 2 (l - z) m (1 + lnz) 

l-u l-u 2{l-u)z l-u 2 (l-u) l-u ' 1 ' 

7vr 2 31n[uz] ln 2 [z u)] / 3 , . . \ , , 

^4-4- — ±-L + —L-Jl +in ^--ln( fi2f )J , (172) 

JA6 = _™ _ >* + *- *Mfi*) , (173) 

7/ ^ -2 

T <^ if 2 - lnii 

■7A7=- — j—, 174 

2 it 2 it 

77T 2 1 

JB4 = i?r ln[wz] + — + -(3- ln 2 [uz]) , (175) 

JB6= 1 -^Hu Z ) + 1 ^°f [UZ] + 7 -^) , (176) 

^1 = ~ . (177) 
2m 

J_B2 = — (—2 in In m — 3 In m + ln 2 u+2 In u In z) , (178) 
2it 

„ 1 77r 2 ln 2 u ln 2 [itz| lnz 3 , , / 3 — 2u lnz 

J£3 = --H — + — lnu +lnu 

\ 'Zuu uu J 

(179) 



u I2u 2uu 2u 2u 2u \ 2uu uu 

(it — 2 \nu — 2u ln[uz]) 
2uu 



JE6 = — !— Uu + (3 - u) hxu - ln 2 u - u lnz - 2 In ft lnz + in (u + 2 Inn)) , (180) 
2uu 

n 2 2 (2-3z)lnz l n fz\ 
JFQ=— + ^- y - —I + -S(-J , (181) 

24u u Zuz u V z / 

J HI = — (4 + in - ln(u z)) , (182) 
2 u 

Jm = + i!LM _ , (i 8 3) 

12 u 2u m 
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JH6 = -L ( 1 + m - ln(u z) ) , (184) 



JCl = ReJCl + mlmJCl, (185) 

(186) 



z (1 — 3 u — z) u 2 z , " // 



j2 ^ 

Im JCl = — tt 5- In 

2(1 -u-z) 2 (1-u-z) 



z J 



u z l(u, z) lnu I u u(l + u — z) 



ReJCl = - — 1 + 



(1-u-z) 2 2 \ u 2 z (l-u-z) 
1 u 2 u 



2 (l-u-z) 2 2(1- u-z) 
1 u u \ (1 — uz) 



2uz (1 — u — z) 2 2u(l — u — z) 



T ™ 7T 2 / 4 n \ 1 2-uV /3z-l., 

JC2 = — 3 - - - — 5- lnu + — lnz- 

21\uz I U 2u 2 z / \ 2uz ' 



In z+ 

ln(l-uz) , (187) 



3 2-u 2u-l 



2u 2u 2 z 2uzu 2 



ln(l — uz)+ 



1 - -4= ] Li 2 («) + ( 1 - 4 J Li 2 (z) - ( 1 - -L ) Li 2 (zu)+ 

HMZ / \ UZ J \ UUZ ' 



u/ \ uz / \ \z/ \ 1 — zm 



JC3 = ReJC3 + i7rImJC3, (189) 

ImJC3=^ 2+lT7^ r -\n(uz) , (190) 

2 (1-u-z) 2 2(1 -u-z) (1-u-z) 3 \z) y >' K 1 



t. T ™ , 7tt 2 u 3 I(u, z) / 3 u 2 u \ , 
RcJC"3 = 1 o - ^ o + ^ r lnw+ 



12 (1-u-z) 2 ^2 (1-u-z) 2 2{1- u-z) 

/ „ 3 u 2 (z- l + 3u)\ lnz 1, 2r 1 
+ 3 - 4u - - + — i > — + - In 2 [uz] + 



3u 3 - 4u (2 - u)u 2 2u(2u - 1) - u 3 
2u 2uzu 2 u(l-u — z) 2 2u 2 (l — u-z) 



(1 — uz) / /z\ / uz 



uz 



5 



1 /2-u 1 3-2u\ _ / 3 1 , 
JC4 = - + — o ^ In z - — 5- (1 - zu) In 1 - zu 

2 V uz 2uz 2 2u J \2uz 2uz 2 u} K 

(1 - zu) / /z\ / uz 



uz \ Vz/ \1 — uz 
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JD1 = RcJDl + inlmJDl, (193) 

lmJDl = — + Z ——\n(-\ , (194) 

u-z (u-zf \u) 

ReJDl = -^-^l+(^-±) ln ^ + i^L_ (l-^)Ml-^) _l s m (195) 
24m u — z \u— z uz J u — z uz(u — z) u \uJ 

TT ^„ 7r 2 /z — 5 — 3uz\ lnu lnz (1 — uz) , . . 
JL>2 = — + ^ '- In 1 - uz] + 

24u \ z J uz zu uzuz 

1 + UZ . /U\ \ — Uz( (Z\ / UZ 

— Li 2 u + Li 2 z - Li 2 uz - S — S - - S 

zu \u) zu \ \zl 



JD3 = RcJD3 + iTTlmJD3, (196) 
ImJD3 = -L + — !— + (4 + 7 ro + — ^— J lnu+ f 2 - - - — j lnz (197) 



2u u — z \u (u — z) 2 u—z J \ u (u — z) 

9 /7 1 \ / 1 — uu uz u \ N / 1 1 

7T 2 — + — 



RcJD'i = + + - - + _ J I(u, z") - ^- + — - ] h, „+ 

1 ~ zu , n t f u 1 1 \ t - 1 , 2 r— 1 

In [1 — zu] — I — + — + ^— I in z — — In [zuj — 



uz(u — z) \2u (u — z) uz 

{z — uz) ( z\ | (1 — uz) Q f uz 
- o 



m+^^S , (198) 

uz Vz/ uz \ 1 — uz / 



T ^ 5 z,_ 1 — uz//z\ / uz 
J£>4 = - H — In z H — I S ( — ) — S 

2 z zu V \zJ 



1 — uz 



JD6 = ReJL>6 + inlmJD6 , (199) 

z z (1 — u (u + z)) , /u\ ,„„„s 

Im JD6 = — + — - — In - (200 

u{u-z) u(u-z) 2 v«y 

ReJ£>6=— +(-- — + -^—\l(u,z)- + )\n(uz) + (1 ~ UZ \ ln(l - uz), (201) 
2u \u u u — z/ \u u—z) uz\u — z) 



where for brevity we used new notation S(x) : 

S(x) = 1 In 2 (1 + x) + [ — ln(l + xa) = Li 2 f -^- > ) + In 2 (1 + x) . 
2 Jo a \l + x J 



All other functions have been defined in the text. 



(202) 
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